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Abstract 

A simple derivation of the static Gross-Pitaevskii (GP) equation is given 
from an energy variational principle. The result is then generalized heuristi- 
cally to the time-dependent GP form. With this as background, a number of 
different experimental areas explored very recently are reviewed, in each case 
contact being established between the measurements and the predictions of the 
GP equations. The various limitations of these equations as used on dilute inho- 
mogeneous condensed Boson atomic gases are then summarized, reference also 
being made to the fact that there is no many-body wave function underlying the 
GP formulation. This then leads into a discussion of a recently proposed inte- 
gral equation, derived by taking the Bogoliubov-de Gennes equation as starting 
point. Some limitations of the static GP differential equation are thereby re- 
moved, though it is a matter of further study to determine whether a correlated 
wave function exists as underpinning for the integral equation formulation. 

1 Introduction 

The purpose of this Chapter is to review recent progress in both theory and exper- 
iment in the area of inhomogeneous condensed bosons. Both statics and dynamics 
will be referred to in the course of this review. However, the field cited above is 
now vast, and therefore we shall select, in both experiment and theory, the areas in 
which our own contribution lie. 

With the above brief outline, we next emphasize that, since the Gross-Pitaevskii 
(GP) equations, to be given explicitly in Section [21 below, have been very valuable in 
providing a theoretical framework, albeit quite approximate, for the interpretation of 
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an extensive body of experimental data, we shall take the GP equations as a focus for 
the present Chapter. However, to avoid repetition, we cite here two major reviews 
which should be consulted by the reader who requires more technical background 
to the basic theoretical arguments, and application techniques, underlying the GP 
equations. 



The first of these is bv lDalfovo et ai.l(|l999h . This provides an authorative review 
of the phenomenon of Bose-Einstein condensation (BEC) of di lute gases in h armonic 
traps from a theoretical perspective. For the ideal Bose gas. lDalfovo et all treat (i) 
finite size effects, (ii) the role of dimensionality, and (Hi) the thermodynamic limit 
of trapped Bosons at finite temperature. In addition, again for the ideal Bose gas, 
nonharmonic traps and adiabatic transformations are also treated. So in the present 
Chapter we shall take their treatment of the ideal Bose gas in traps as assumed 
background. 

The second of the reviews referred to above is of very different character, but 
again provides invaluable background for this present Chapter. Thus the focus of 



the review by IMinguzzi et all ((200 4:) is on numerical methods in use for atomic 
quantum gases with app lication to Bose -Einstein condensates, our dominant focus 
in this Chapter, but also iMinguzzi et all cover ultracold Fermions. Thus the reader 
interested specifically in numerical techniques used to solve practical problems in 
inhomogeneous assemblies of condensed Bosons should consult the extensive treat- 



ment of IMinguzzi et al. 



( 2004h . 



With this background, the outline of the present Chapter is as follows. In Sec- 
tion |21 both static and dynamic GP equations are set up by what we consider to be 
the simplest physical considerations. One example to illustrate the time-dependent 
GP Eq. © below, which will be su mmarized also in Section |2l is drop emission 
from an optical lattice under gravity ( Anderson and Kasevichl . 19981 : Chiofalo et al. , 
1999), while a somewhat generalized static GP equation is utilized in Fig. ^ which 
shows a vortex array for a rotating Bose-Einstein condensate. These two exam- 
ples lead into Section 03 in which selected experiments are both described and also 
brought into contact with the GP equations. In particular, we shall focus on opti- 
cal lattices, atom chips and magnetic microtr aps, and the real i zation of Josephson 



Junction arrays of Bose-Einstein condensates ([Cataliotti et all . l200ll ) . We will also 



briefly survey the dynamics of a BEC expanding in a moving ID optical lattice, 
and some experiments beyond the GP equations. Se ction H] is then a very brief 
discussion, following the arguments of Leggettl (|200J ). of some limitations of the 
GP equations exposed by first-principles arguments. This leads into Sections El and 
H3 in which an integral equation trans c endin g the static GP Eq. © is set up, fol- 



lowing the proposal bv lAngilella et all (J2004). The Chapter concludes with a brief 
summary and some selected proposals for directions in which further studies should 
prove fruitful (Section EJ. 
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2 Gross-Pitaevskii (GP) equations 



The equat ions which p r ovide a focus f or the pres e nt rev iew were proposed indepen- 
dently by iGrossI (|l96ll . 1963) and bv IPitaevskiil l|l96ll ). However, it i s fair to say 
that t he mean- field description of a dilute Bose gas goes back at least to lBogoliubov 
(1947). The central point in his study was in separating out the condensate contri- 
bution. The generalization of the Bogoliubov approach to embrace the situations 
where there is both inhomogeneity and time-dependence is the aim of this section 
of the present review. 



2.1 Static GP equation 

The natural starting point then is to treat first the static, or time-independent 
case, for inhomogeneous condensed Bosons. This will lead t o the so-called static 



GP equation. The a pproach we shall adopt below (see also Dalfovo et al. . 19991 : 



Mm^zi ,t A H is to start from variational minimization, for a Bose-cond~d 
gas in a three-dimensional trap at zero temperature, of the energy functional E[<&\ 
given by 



dr 



2m 



V$| 2 + F ext (r 



|$f + -s|$| 4 



(1) 



In this energy functional, |<I>(r)| is the inhomogeneous density profile, the so-called 
'condensate wave function', $(r) playing the role of the order parameter of the 
Boson assembly. The first term on the right-hand-side (RHS) of Eq. (^Q) represents 
the kinetic energy of the condensate while the second is the potential energy due 
to the external confining potential T4 x t(r). The final term on the RHS of Eq. (^Q) 
requires rather more discussion. It is appropriate in a dilute and cold g 
representation of the atom-atom interaction, say in general U(r — r'), since in this 
case only binary collisions at low energy prove to have relevance. These collisions 
can be characterized, in fact, by a single parameter a, which denot es the s-wave 
scatte ring length, independently of the finer details of U(r — r'). As lDalfovo et al 



(1999) stress, this permits one to adopt an effective interaction given by 



g5{r 



(2) 



(see also Parkins and Walls) . ll99Sl : iLeggettl . l200ll ) . Here the coupling constant g is 
given in terms of the s-wave scattering length a in three dimensions and atomic 
mass m by 

Airh 2 a 

9 = • (3 

m 

It has been emphasized that g is a parameter depe nding on particle density 
( Pieri and Strinatil . [2000) as well as on dimensionality (jChernv and BrandL [2004, 
an d references therein). In particular, its dependence on density has been discussed 
by Pieri and Strinatil ( 2000l ) in connection with the crossover from wek-coupling BCS 
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super conducti vity to strong-co upling; BE condensation ( Nozieres and Schmitt- Rink , 
198,4 see also iRanderial Jl995) for a review). 

Returning then to the functional Eq. it is to be noted that for repulsive 
interactions corresponding to g > 0, the functional is convex and the minimum 
corresponds to the stable ground state. It should be emphasized, for the other case 
when g < 0, that the ground state exists only at weak coupling for a limited number 
of trapped Bosons, so long as the zero-point energy is able to balance the effect of 
attractions and prevent collapse. 

With this short introduction to the assumed energy functional E[<&] set out 
in Eq. (fljl. it is then a straightforward matter to perform the minimization with 
respect to the order parameter $(r). One is then led to the (of course approximate) 
static GP differential equation having the form of a non-linear Schrodinger equation, 
namely 

h 2 



—V 2 + V ext (r)+g\<i>(r) 



2m 



$(r) =/i$(r), 



(4) 



a condensate at equilibrium being at chemical potential [i. 

While we shall consider rather fully below a variety of recent experimental data 
obtained by some of the present authors, let us briefly summarize at this point just 
one app lication of the sta t ic GP equation, Eq. (@J). The example below follows the 
work of iGastin and Duml |l999) who treated a vortex array. This they obtained 
by numerical solution of the two-dimensional differential equation for a rotating 
Bose-Einstein condensate (BEC) under tight vertical confinement. In this rotating 
condensate, described in the rotating frame by an appropriate static GP equation, 
one needs to incorporate an inertial term. Then Eq. can be extended to the form 



7^V 2 + y cxt (r) + <7|$(r)| 2 + ftL, 
2m 



/i$(r), 



(5) 



Q denoting the rotation frequency and L z the an gular momentum c ompo nent along 
the rotation axis. As emphasized in the study of Castin and Duml ( 19991 ). the term 
£IL Z in Eq. © is not diagonal in position or in momentum space, and this implies 
that numerical solution of Eq. © requires special care. 

Starting with different trial states, a single-vortex solution, as well as those 
describing multi-vortex configurations, can be exhibited. In this latter form of so- 
lution, the vortices display a tendency to arrange th emselves into a tr i angul ar ge- 
ometry, as shown in Fig. [T] which is redrawn from Castin and Duml (^)99). As 
stressed by iMinguzzi et al a significant result which has emerged from 

solving the static GP Eq. (J5J) in a cigar-shaped trap i s that a vortex line can 



be bent, as shown in Fig. 17 of Minguzzi et al. ((20041 ) , following the study of 



Garcia-Ripoll and Perez- Garcia! (|200ll ). This finding is in general agreement with 
experimental observations. 
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Figure 1: Shows vortex array as exhibited by numerical solution of the static 
Gross-Pitaevskii Eq. ffl for a rotating Bose-Einstein condensate. Redrawn from 
Castin and Duml (|l999h . 



2.2 Inhomogeneous Bose superfluids: 
dependent GP equation 



dynamics using the time- 



Let us proceed to treat in a somewhat parallel manner the dynamics of inhomo- 
geneous trapped superfluid Boson gases. One can make a heuristic generalization 
of the non-linear Schrodinger equation having the form of the static GP Eq. to 
describe the now time-dependent condensate wave function $(r,i), to read 



■^-V 2 + y cxt (r) + <7|<D(r,t)| ; 
2m 



$(r,t) = ift|U(r,i) 



(6) 



which, of course, reduces to the static Eq. 0) o n mak i ng th e sub stitution $(r,t) = 
$(r) exp(-iu.t/h). Eq. ©, also given by brossl dl96ll . Il96.^ and IPitaevskiil <ll96lh . 
will be discussed in some detail below, in order to gain further insight into the 
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realm of validity. Following Dalfovo et al. ( 1999h . let us write down the many- 
body Hamiltonian H, in second quantization, which describes N interacting Bosons 
confined again by an external potential V ex f This reads: 



H 



dr^(i 



2m 



V 2 + Fext(r) 



9(r) 



+ 



\jj drdr'¥(Y)¥(Y')U(Y - r')*(r')£(r) 



(7) 



where U(r — r') is the two-body interatomic potential already referred to. In Eq. (J7J), 
^(r) and ^(r) are the Bosonic field operators in the Schrodinger representation 
(t = 0, say). Then, it follows that the evolution equation for the operator \f(r, t) 
takes the form 



2m 



V 2 + I4xt(r)+ / dr'^(r',t)U{r-r')^(r',t) 



*(r,<). (8) 



For a very dilute and fully Bose-Einstein condensed gas at T = 0, it is assumed that 
the field operator can be replaced by the classical field <3?(r, t). Taking again the 
description of U{r — r') in Eq. (J2J) with g related to the s-wave scattering length a 
by Eq. ©, one is led back to the time-dependent GP Eq. ©. 

It is known that the time-dependent GP equation, because it omits dissipation, 
is not appropriate to deal with the dynamics of a confined Bose gas when quantum 
depletion or thermal excitations play a significant role. As examples of such sit- 
uations, one may cite condensate formation and decay, phase coherence, damping 
of collective motions and excitations from a non-mean-field ground state. Nev- 
ertheless, the time-dependent GP Eq. © has been found valuable in treating a 
variety of dynamical processes in condensate clouds. One may cite, as examples, 
collective- mode frequencies in harmonic or optical-lattice confinement, interference 
phenomena, dynamics of vortices, propagation of solitons, shock-wave dynamics, 
four-wave mixing, atom-laser output, as well as expansion of a rotating condensate. 

While numerous examples of the use of both GP Eqs. Q and © will be given 
below in Sec. El one imme diate illustration to s how the usefulness of Eq. © will 
be taken from the work of lChiofalo et al. Their extensive numerical study 

of drop emi ssion from an optical la t tice u nder gravity was motivated by the ex- 
periment of lAnderson and Kasevichl (jl998h . In their experiment, an almost pure 
BEC was poured from a magneto-optic trap into a vertical optical lattice which was 
produced by a detuned standing wave of light from the counter-propagating laser 
beams. The gravitational force tilts the lattice potential and drives tunnelling from 
well states to the continuum. Interference taking place between coherent blobs of 
condensate at different lattice sites leads to the appearance of falling drops, which 
can be regarded as coherent matter-wave pulses by analogy with a mode-locked pho- 
ton laser. Analogous to the laser cavity is the Brillouin zone in momentum space, 
so that the modulation interval for the pulsed emission of drops equals the period of 
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Figure 2: Plots shown relate to condensate drop emission from an optical lattice 
under gravity. Actually shown is the density profile of the condensate drops after 
4.6 ms as a function of a (suitably scaled) distance z in the non-interacting limit. 
The central condensate around z = has been substracted out. Different curves 
plotted are characterized by increasing well depths, going from shallow ( interme diate 



thickness curve) to deep (thickest curve). Redrawn from Chiofalo et al. (1999). 



Bloch oscillations. The time interval between successive drops was 1.1 ms in both 
the experiment and the numerical study, which is the period of Bloch oscillation 
for a quasi-particle (the whole coherent condensate) driven by the constant force of 
gravity through a periodic array of potential wells having a lattice period of one-half 
of the laser wavelength A. 

Primary attention was paid bv lChiofalo et al\ (|l999h in their numerical investiga- 
tions to the part played by atomic interactions in this example of coherent transport. 
In agreement with band-structure theory of transport in periodic structures, the pe- 
riod of drop emission depends only on the intensity of the drive and the lattice 
spacing. The size and shape of the emitted drops are instead determined by the 
magnitude of the lat tice barrier and the interatomic forces. In Fig. [2 redrawn from 



Chio falo et al\ ((199 9 ) , the density profile of the drops emitted from a non- interacting 
BEC after 4.6 ms is depicted as a function of a suitably scaled distance z. The main 
condensate around z = has been subtracted, the curves being characterized by 
values of the lattice barrier height. In the experiment, the latter is determined by 
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the intensity of the optical beams. 

When the constant drive is supplemented by a monochromatic oscillating drive, 
it turns out that the equations governing this form of coherent BEC transport in the 
linear non-dissipative regime can be mapped onto those for the ac superconductin g 
current flowing across a weak- link Josephson junction ( Chiofalo and TosH l200lh . 
The voltage drop across the junction is, in essence, replaced by the product of 
the constant component of the force times the lattice spacing, and the frequency 
of the Bloch os cillations is in reson ance with integer multiples of the oscillating 
drive frequency. iBurger et al. ( 200 il l have carried out the relevant experiment on a 
BEC confined in a magnetic trap plus an optical lattice, the oscillating force being 
generated by a rapid shift of the centre of the trap. 

With this introduction to both static and time-dependent GP differential equa- 
tions given in Eqs. @ and @ respectively, we turn next to discuss a variety of 
recent experiments by some of the present authors, and the relevance of these GP 
equations to their interpretation. 



3 Selected experiments brought into contact with GP 
equations 

The first experimental realizatio n of a Bose-Einstein condensate (BEC) in a di- 
lute gas of rubidium atoms (see llnguscio et al. , ll 99fiL as a general reference) has 
marked the birth of a completely new field of optics, that of coherent matter wave 
beams. Indeed, it was immediately recognized that a BEC is in very many aspects 
the analogue of an optical laser for atoms. After the first pioneering experiment at 
MIT where a pulse d atom laser was cr eated by simply switching off the confining 
magnetic potential ( Mewes et al. . 199?t ). other more refined experiments followed. 
The first continuous wave atom laser was realized i n Munich by coupl ing atoms 



out of a magnetic trap using a radio frequency field ( Bloch et al. . 1999). while in 



Yale a BEC was loa ded into a vertical standing w ave to produce the analogue of 
a mode- locked laser ( Anderson and Kasevichl . 199£ ). Coherent matter w ave beams 



were soon shown to produce interference fringes bv I Andrews et al, 
ble slit and mult iple beam interferometers were soon realized i 



dl997h. and dou 



Fort et al. 2001 



Pedri et all 12001). On the other hand atoms, at variance with photons, interact 



even in the dilute gas limit. The nonlinearity introduced by atomic interactions in 
coherent matter waves propagation through vacuum was readily shown to be the 
same of the third order susceptibility for electromagnetic waves propagating through 
a nonlinear crystal. It was soo n demonstrated that i t was indeed possible to ob- 



serve matter wave amplifi c ation ( Kozuma et all 19991 : Inouve et al. , Il999l ) and four 



wave mixing (|Deng et ail Il999h . More recently, the intrinsic nonlinearity of con- 
densates was successfu lly exploited in the o bservation of dark ( Burge r et al. , 1999; 



Denschlag eTaA. l200(t) and bright solitons (|Khavkovich et al. 



2002). The last re- 



sult could also be obtained by controlling the dispersion of matter waves using a 
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periodic potential ( Eiermann et al . 2003h 



3.1 Standard experimental procedures 

A dilute gas of particles can be cooled to a temperature such that the de Broglie 
wavelength associated to each particle becomes larger than the mean interparticle 
distance. At these temperatures the quantum statistics of the particles fully domi- 
nates the behaviour of the gas. In a trapped atomic gas of bosons all the particles 
will tend to occupy the trap state with the largest population in a very similar 
way to photons in a laser cavity being pulled to the mode with the highest gain. 
In thermal equilibrium the state with the highest occupancy is the ground state 
of the trap. In a gas of interacting bosons, however, the ground state of the sys- 
tem will not necessarily be the ground state of the potential holding the atoms. 
In a dilute ultracold atomic gas of bosons only binary collisions are allowed and 
the system is conveniently described by the Gross-Pitaevskii equation @ for the 
atomic density amplitude. In particular, it is possible to change the value of the 
s-wave scattering length a by different orders of magnitude and eve n change i ts sig n 
by tuning a magnetic field thanks to Fano-Feshbach resonances ( Ingusciol . Ef 



The parameter a strongly influences the properties of the condensate. Condensates 
with negative a are unstable but can form bright solitons. On the other hand, the 
ground state of condensates with repulsive interactions (a > 0) significantly differs 
from the ground state of the confining potential. Indeed, for most experimentally 
realized condensates the nonlinearity totally dominates dispersion, i.e the interac- 
tion energy is much larger than the kinetic energy. The condensate wavefunction 
then takes the so called Thomas-Fermi shape which has the same symmetry of the 
trapping potential an d, for a condensate in a harmonic trap, is an inverted parabola 



I Dalfovo et all \U 

The standard experimental procedure for the creation of a Bose-Einstein con- 
densate in a dilute atomic gas starts with laser cooling of an atomic vapor in 
ultra- high vacuum con ditions. This first step, performed in a magneto-optical 
trap (R aab et all Il987l ). takes the atoms from a phase space density of 10~ 20 at 



room temperature to 10 5 below 100 fiK. At this phase space, density l aser cooling 
stops essentially because of spontaneou s light scattering from the atoms 



ser c ooling 



Cohen-TannoudiiL Il998l : IPhillipsl 1998). For this reason, the following cooling step 



has to be performed in a non-dissipative trap created either by a magnetic field or 
by a very far-off resonance laser beam. In non-dissipative traps, cooling is achieved 
by removing the high energy tail of the atomic distribution and by letting the atoms 
thermalize via binary collisions. Removal of atoms is realized either by reducing 
the trap depth in optical traps or by RF-induced transitions to untrapped Zeeman 
sublevels in magnetic traps. 

Let us concentrate on magnetic traps. The interaction energy of an atom in 
a magnetic field is E = — /x • B, where n is the atomic dipole moment. Of the 
atomic ground sublevels, only low-field seeking states can be trapped in a magnetic 
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field minimum. If the atomic motion is not very fast, the atomic dipole adiabatically 
follows the magnetic field. Therefore the energy is proportional to the modulus of the 
field. However, if the field variation is very rapid, adiabaticity cannot be maintained 
and atoms are lost from the trap. To avoid these so called Majorana spin- flips, it is 
necessary to avoid a zero of the magnetic field in the trap. Conventional magnetic 
traps are realized with a few centimeter size coils carrying more than 100 Ampere 
of current. In this way traps with 10 — 100 Hz oscillation frequencies are realized. 

The cooling procedure, known as evaporative cooling, strongly relies on atomic 
collisions and can be very different for different atomic species. In particular, 
since the most probable collisions at temperatures below 1 mK are sphericall y 



symmetric binary collisions, due to the Pauli principle ((DeMarco et all llf 
this method cannot work for spin-polarized fermions as those trapped in a mag- 
netic t rap. Indeed, fermions were cooled to degeneracy either using spin mix- 



tures (foeMarco and Jirl 11 999) or v i a col li sions with a bos o nic a t omic species 



dTruscott et all l200ll: ISchreck ' et all US]; iHadzihahic etUl El Iftmti al 



200 



r 



3.2 Optical dipole potentials 

Matter wave beams can be manipulated in very much the same way as optical 
beams. However, the role of matter and electromagnetic fields is totally reversed in 
the field of atom-optics. Indeed, condensates can be manipulated with conservative 
potentials created either by far-off resonance laser beams or by magnetic fields. In 
the first case the potential is obtained via the interaction of the induced atomic 
dipole with the electric field of the laser. This dipole potential is dependent on the 
laser intensity and detuning and for a two level atom in interaction with a far-off 
resonance beam can be written as 



3ttc 2 r 

4 



where loq is the atomic resonance frequency, V the natural linewidth of the atomic 
transition, A the laser detuning from resonance, and I(r) the intensity of the laser- 
beam. From Eq. Q we note that when the laser detuning is negative, atoms are 
pulled towards the region with the highest laser intensity. On the other hand, when 
the detuning is positive, atoms are expelled from high intensity regions. This can 
be used to create very different kind of potentials with a single tunable laser beam. 
For example, it is possible to create atomic waveguides either by using collimated 
gaussian laser bea ms with nega t ive det uning, or collimated hollow laser beams with 
positive detuning ( Bongs et al. . 2001b! ) . 



With a sheet of light created by rapidly moving a collimated laser beam with 
an acousto-optic modulator, the group of IQO-Hannover has been able, by simply 
varying the laser beam intensity, to create an atom mirror, a beam-splitter, or a 
phase-shifter (see Fig. 0J. The laser detuning for these experiments was positive. 
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User power 



Figure 3: Scheme of the IQO experiment. The condensate falls under the effect of 
gravity and bounces off the potential generated by the laser light. 



The experimen t was performed by c reating a condensate via evaporative cooling in a 
magnetic trap ( Inguscio et al. , 19981 ) and then dropping it under the effect of gravity 



onto the sheet of light. When the laser intensity was very high, the atoms did not 
acquire sufficient kinet ic energy dur i ng th eir fall to go over the dipole potential and 



were totally reflected ((Bongs et all [1999). When the laser intensity was reduced, 
part of the atoms could go over the dipole potential and the system realized a beam- 
splitter. By further reducing the laser intensity all the atoms were able to go over 
the potential but were retarded w ith respect to free fall in analogy to an optical 
phase-shifter iBongs et all 12001a), 

3.3 Atom chips and magnetic microtraps 

These examples demonstrate that laser radiation is a very versatile tool, since it 
is possible to create many different potentials simply from the interference of laser 
beams coming from different directions as will be discussed below. On the other 
hand, it is also possible to manipulate atoms with magnetic field gradients. Magnetic 
traps are, for neutral atoms, a particularly versatile class of these manipulation 
methods as they can be used for any atomic species with a magnetic moment and 
they can produce conservative potentials also for very long times. Techniques to trap 
and manipulate atoms with magnetic fields once integrated with surface deposition 
techniques, either lith ographic or of other kind, realize what is termed an atom chip 
(Fol man et al. Atom chips are based on the possibility of creating a 2D- 



quadrupole magnetic field close to a current carrying wire by compensating the field 
generated by the wire B = at the height zq with a constant magnetic field. If 



12 



G. G. N. Angilella, S. Bartalini, F. S. Cataliotti et al. 




* + 




wire field 



bias field 



Figure 4: Using a bias field to compensate at a point zq the magnetic field generated 
by a current carrying wire, it is possible to create a 2D-quadrupole trap. 

the magnetic moment fi of the atom remains aligned with the magnetic field the 
resulting potential for the atoms can be approximated as 



where we have assumed the wire to be along the x axis. With a current of just 0.4 A 
in the wire and a constant field of 35 G, it is possible to create a waveguide with 
a confining frequency of 10 kHz for 87 Rb much larger than conventional magnetic 
traps. 

The wire guide represents the building brick for magnetic microtraps. Indeed, 
by bending the wire in a U shape as illustrated in Fig. 0(b), it is possible to create 
a 3D quadrupole. As already pointed out, this trap configuration have high atomic 
losses due to Majorana spin-flips. In order to avoid these losses, it is better to create 
a harmonic trap with a filed minimum different from zero as can be done by bending 
the wire in a Z shape [Fig. (a)]. The bias field can also be realized in a planar 
configuration using three parallel wires with opposite currents as in Fig. 0(d). This 
configuration can then be used to create more complicated structures as in Fig. 0(c), 
where a possible analogue of a SQUID is illustrated. 

These devices represent one of the most promising schemes for coherent atom 
optics and may be the basis for a totally new class of integrated sensors and quan- 
tum logic instruments. Indeed, thanks to component miniaturization, in atom chips 
it is possible to reach huge field gradients (above 1 Tesla/cm) with 1 Ampere cur- 
rents, of the same order of magnitude as those used in electronic circuits and a few 
orders of magnitude lower than those employed in conventional apparatuses. Fur- 
thermore, the substrates employed are compatible with the ultra-high vacuum tech- 
nology needed for atom cooling. Nowadays the microstructures used experimentally 
are rather simple and the conductors sizes are of the order of 10 /mi with a total chip 
size of a few cm 2 . However, they are capable of substituting experimental systems 
currently spread in about 1 m 3 . Integrating many elements to control atoms onto a 
single device, an atom chip, will make atom physics experiments much more robust 
and simple. This may allow much more complicated tasks in atom manipulation to 
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Figure 5: Various microtrap configurations: (a) Z-Trap; (b) [/-Trap. Both the Z- 
Trap and [/-Trap need an external bias field to achieve a 3D-confining potential for 
cold atoms; (c) A ring guide with two diametral potential barriers, working with 
cold atoms like a SQUID works with superconductors; (d) A linear guide completely 
realized in planar geometry, without any bias field. 



be performed in a way similar to how integration of electronic elements has allowed 
the development of new powerful electronic devices. Ideally, one desires routine sin- 
gle atom state selective loading, preparation and manipulation. All will be achieved 
with minimal heat load and power consumption. Potentials with sizes smaller than 
the particle de Broglie wave-length will allow tight traps with large energy level 
spacing. The large level spacing reduces the probability of the environmental noise 
to induce unwanted excitations. Consequently, coherent manipulation will be more 
stable. Using well developed nanofabrication technology from microelectronics to 
build the atom optics will allow integration of many atom optical devices into com- 
plex quantum networks combining the best of two worlds: the ability to use cold 
atoms, a well controllable quantum system, and the immense technological capabil- 
ities of nanofabrication and microelectronics to manipulate the atoms. Experiments 
based on atom chips have allowed the demonstration of trapped a t om int erferom- 
eters ((Hansel et aZl l200lJ>. of atomic co nveyor belts ((Hansel et ail l2001bh . and of 



atomic waveguides (jFolman et all l200 Ct). Currently one of the main concerns in 
atom chips is the possibility of maintaining a coh erent sample close to a surface . 
After the achievement of BEC with these devices ( Hansel et ai . 200 lal : Ott et al. 
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200lh . it wa s discovered a f ragme ntation of the atomic cloud when brought close to 



"surface JFortagh et a/I EoO^. This fragmen tation was ver y recently attributed 



to the thickness fluctuations of the conductors ( Esteve et al. . 20041 ). which causes 



the current to deviate from a straight line producing unwanted a xial magnetic fields. 



Other experiments are concentratin g on surface induced losses ((.Tones et all |2 



Harber et all 1200.4 iLin et all l2004h . however it was possible to realize an atomic 



clock starting from a BEC in an atom chip ( Trcutlei n et all |2 



3.4 BEC in periodic potentials 

We now describe some experiments performed on BECs in periodic potentials, 
that well demonstrate the flexibili ty of optic a l pot entials in the manipulation of 
condensates. In their experiment, iFort et al. (2000) load cold 87 Rb atoms in the 
\F = Imp = —1} state from a double magneto-optical trap system into a Ioffe 
type harmonic magnetic trap. The trap is cylindrically symmetric wit h an axial 



freque ncy of uj x /2tt = 9 Hz and a radial frequency of ui±/2ir = 92 Hz. IFort et al 
(2000) ;hen cool the atoms in this trap via rf-forced evaporation until they reach 
a temperature below 300 nK, just above the critical temperature for condensation, 
which is around 150 nK dependent on the number of atoms we load in the trap. At 
this stage of evaporation they suddenly switch on an optical standing wave formed 
by retroreflecting light from a laser blue detuned of ~ 3 nm with respect to the D\ 
transition at A = 795 nm. The laser beam is aligned horizontally along the axis of 
the magnetic trap and, since the laser beam waist is much larger than the atomic 
cloud transverse size and does not produce any appreciable radial force, forms an 
array of disk shaped traps together with the magnetic potential. 
The potential is therefore 



V 



-m(uj x x + to±_r ± ) + sEr cos (2irx/\). 



(11) 



The optical potential is given in units s of the energy Er = h 2 /2mX 2 gained by 
an atom (of mass m) absorbing one lattice photon corresponding in rubidium to 
a temperature of ~ 170 nK. In these experiments the optical potential could be 
varied up to s = 15. When s > 1 the atoms are confined in an array of classically 
independent traps since the optical potential barriers are much higher than the 
thermal energy of the atoms. Indeed it is possible to "freeze" the degree of freedom 
associat ed with the motion along the axis and study condensation in this quasi 2D 
system (B urger et al. . 2002). After switching on the laser light we continue the 
evaporation ramp until the desired temperature is reached. This ensures that the 
atoms reach the equilibrium state in the combined trap. When we evaporate to 
well below the critical temperature, so that no thermal fraction is experimentally 
visible, we typically obtain ~ 200 condensates separated by a distance of A/2, each 
containing ~ 1000 atoms. Due to the blue detuning of the laser beam the atoms are 
trapped in the nodes of the standing wave reducing light scattering below ~ 1 Hz. 
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When the height of the optical barriers is much larger than the condensates 
chemical potential we are justified in describing the condensate as a sum of wave- 
functions localized in each potential well: 

*o(r)= exp[-(z-£;A/2)V2a 2 + i0 fc ]^^Yl- 7 ^y (12) 

k=o,±i...±fc M 9 V y K± )k/ 



where (R±)k = \l'^P t -k/' muJ '±_ ls the radial size of the k-ih condensate, g depends on 
the scattering length a through Eq. ©, while \ik = ^mu> 2 d 2 (k 2 M — k 2 ) plays the 
role of an effective /c-dependent chemical potential. The value of ku is fixed by the 
normalization condition N = to 



^ 2 27iu; / 15 a d 
mu) 2 d 2 vSy 7 ^" a,ho & 



with uj = (uJxUJ 2 ^) 1 / 3 the geometrical average of the magnetic frequencies, a^ G = 
\Jhjmijj is the corresponding oscillator length. From the above equations one also 
obtains the result N k = N (l - k 2 /k 2 M ) 2 with N = ^N/k M . Equation (JEJ) gener- 
alizes the well known T h omas -Fermi results holding for magnetically trapped con- 
densates (Dalfovo et al. , 1999) to include the effects of the optical lattice. This 



generalization is justified by the fact that the optical confinement along the optical 
lattice is much stronger than the magnetic potential therefore it is more suitable 
to use a harmonic approximation for the wave-functio n along the x-dire ction and a 
Thomas-Fermi approximation in the radial direction (jPedri et al. 



When the atoms are released from the combined trap they spread out and over- 
lap producing an interferogram which will depend on the relative phases (ftk of the 
individual condensates. In Fig. (A), we show a typical image of the cloud taken 
after an expansion time t e xp = 29.5 ms, corresponding to a total number of atoms 
iV ~ 20000 and to an optical potential s = 5. The image shows a clear structure with 
three interference peaks separated by 2h/mXt exp , i.e. by the distance corresponding 
to the reciprocal of the lattice constant. We remark that, differently from the case 
of two separated condensates, interference fringes appear only if the initial configu- 
ration is mutually coherent. In other words, since one single interference experiment 
with an array of condensates is equivalent to averaging a series of interference exper- 
iments with two condensates, an interference pattern will appear only in presence 
of a fixed relative phase between condensates belonging to consecutive wells. What 
is locking the phase difference across the array of BEC is tunnelling through the 
optical barriers, in a classical picture no interference peaks would arise. The width 
of the central peak (n = 0) of the interferogram is of the order Ap x ~ h/mR x t exp 
where R x ~ kud is half of the length of the whole sample in the x-direction. The 
occurrence of these peaks is the analogue of multiple order interference fringes in 
light diffraction. 
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distance (nm) 



Figure 6: Left: (A) Absorption image of the density distribution of the expanded 
array of condensates. (B) Experimental density profile (crosses) obtained from the 
absorption image (A) integrated along the vertical direction. The wings of the cen- 
tral peak result from a small thermal component. The continuous line corresponds 
to the calculated density profile for the expanded array of condensates for the ex- 
perimental parameters (s = 5 and t e xp = 29.5 ms). Right: Experimental (circles) 
and theoretical (triangles) values of the relative population of the n = 1 peak with 
respect to the n = central one as a function of the intensity factor s of the optical 
potential V op t. 



The relative population of the n ^ peaks with respect to the central one 
(n = 0) obeys the simple law 

P n = exp(-16vr 2 nV/A 2 ), (13) 

holding also in the presence of a smooth modulation of the atomic occupation num- 
ber Nk in each well. Equation 1)13(1 shows that, if a is much smaller than A/2, 
the intensity of the lateral peaks will be high, with a consequent important lay- 
ered structure in the density distribution of the expanding cloud. The value of a, 
which characterizes the width of the condensates in each well, is determined, in first 
approximation, by the optical confinement. By using a numerical minimization of 
the energy we can determine the relative population P n of the n = 1 peak as a 
function of the intensity factor s. This is shown in Fig. H3 (right) together with the 
experimental results. The good comparison between experiment and theory reveals 
that the main features of the observed interference patterns are well described by 
this model. 
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3.5 Josephson Junction array with BECs 



In the preceding subsection we have shown that the BECs produced in the combined 
trap are phase-locked by the tunnelling of atoms through the optical barriers. The 
system indeed realizes a one-dimens ional array of Jo s ephso n Junctions (JJs), as we 



wish to demonstrate in this section ( Cataliotti et al. . 2001) 



A Josep hson junction is a sim p le dev i ce made of two coupled macroscopic quan- 
tum fluids ( Barone and Paternol . 19821 : Baronel . 2000h . If the coupling is weak 
enough, an atomic mass current I flows across the two systems, driven by their rel- 
ative phase A(j) with a limiting current I c , the "Josephson critical current", namely 
the maximal current allowed to flow through the junction. The relative phase dy- 
namics, on the other hand, is sensitive to the external and internal forces driving 
the system being driven by the chem ical potential difference between the two quan- 



tum fluids (|Tillev and Til lev. 1990|). The arrays of JJs are made of several simple 
junctions connected in various geometrical configurations. In the last decade, such 
systems have attracted much interest, due to their potential for studying quan- 
tum p hase transitions in systems w here the external parameters can be readily 
tuned ( Fazio and van der Zantl . 1200 ll ) . Recently, the creation of simple quantum- 
logic units and more complex quantum computer schemes have also been discussed 
(jMakhlin et adkOQlh . 

The condensates in two neighbouring sites of the array have a significant in- 
teraction via the tunnelling through the barrier, we can therefore rewrite the time- 
dependent Gross-Pitaevskii equation ©, normally used to describe weakly interact- 
ing conde nsates, as a discrete non-linear Sc hrodinger (DNLS) equation in a parabolic 
potential ( Trombettoni and Smerzl 200lh . We remark that the model used here is 
one dimensional, as we are now only concerned with the motion along the array. 
We write the condensate order parameter as ^(x,t) = Ylj \J Nj (t) e l< ^ 3 W 5^ (x) and 
obtain 



th- 



at 



-K(lp n -i + Vn+l) + (e n + A | I )lpn, 



(14) 



where e n = f2n 2 , Q 



mm: 



2 I A 



(|) , A = g N T fdx$j. The tunnelling rate is K 



j 

We observe that the wavefunctions <&j, as 
well as the tunnelling rate K, depend on the height of the energy barrier. 

In the ground state configuration the Bose-Einstein condensates are distributed 
among the sites at the bottom of the parabolic trap. If we suddenly displace the 
magnetic trap along the lattice axis by a small distance ~ 30 /im (the dimension of 
the array is ~ 100 /xm), the cloud will be out of equilibrium and will start to move. 
As the potential energy that we give to the cloud is still smaller than the inter-well 
barrier, each condensate can move along the lattice only by tunnelling through the 
barriers. A collective motion can only be established at the price of a well definite 
phase coherence among the condensates. In other words, the relative phases among 
all adjacent sites should remain locked together in order to preserve the ordering of 
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the collective motion. The locking of the relative phases will again show up in the 
expanded cloud interferogram. 

For not too large displacements, we observe a coherent collective oscillation of 
the condensates, i.e. we see the three peaks of the interferogram of the expanded 
condensates oscillating in phase thus showing that the quantum mechanical phase 
is maintained over the entire condensate (Fig. [7[ left). In the top part of Figure [7| 
we show the positions of the three peaks as a function of time spent in the com- 
bined trap after the displacement of the magnetic trap, compared with the motion 
of the condensate in the same displaced magnetic trap but in absence of the op- 
tical standing- wave. (We refer to this as to "harmonic" oscillation.) The motion 
performed by the center of mass of the condensate is an undamped oscillation at 
a substantially lower frequency than in the "harmonic" case. We remark that in a 
thermal cloud, although individual atoms are allowed to tunnel through the barri- 
ers, no macroscopic phase is present and no motion of the center of mass should be 
observed. The center of mass positions of the thermal clouds are also reported in 
Fig. (left) together with the "harmonic" oscillation of the same cloud in absence 
of the optical potential. As can be clearly seen, the thermal cloud does not move 
from its original position in presence of the optical lattice. 

The current flowing through the junction between two quantum fluids has a 
maximum value, the critical Josephson current I c , which is directly proportional 
to the tunneling rate K. The existence of such a condition essentially limits the 
maximum velocity at which the condensate can flow through the inter-well barriers 
and therefore reduces the frequency of the oscillations. As a consequence, we ex- 
pect a dependence of the oscillation frequency on the optical potential through the 
tunneling rate. If we rewrite the DNLS equation (|14|) in terms of the canonically 
conjugated variables population/phase and use collective coordinates, we arrive to 
a phase-current relation 



where £(i) is the center of mass of the array and A(j)(t) the relative phase across 
the junction. We remark that, in the regimes we are considering, the current- 
phase dynamics does not depend explicitly on the interatomic interaction. However, 
it is clear that the non-linear interaction is crucial on determining the superfluid 
nature of the coupled condensates, by locking the overall phase coherence against 
perturbations. 

From Eqs.ED we can see that the small amplitude oscillation frequency uj of the 
current I = Nt4z£ gives a direct measurement of the critical Josephson current I c = 
2KNx/h and, therefore, of the atomic tunneling rate of each condensate through the 
barriers. The critical current is related to the frequency u of the atomic oscillations 
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Figure 7: Left: (A) Center of mass positions of the three peaks in the interferogram of the 
expanded condensate as a function of the time spent in the combined trap after displacement 
of the magnetic field. Up and down triangles correspond to the first order peaks, filled circles 
to the central peak. Open circles show the center of mass position of the BEC in absence 
of the optical lattice. The continuous lines are the fits to the data. (B) Center of mass 
positions of the thermal cloud as a function of time spent in the displaced magnetic trap 
with the standing wave turned on (filled circles) and off (open circles). Right: The ratio 
of the frequency of the atomic current in the array of Josephson junctions to the harmonic 
trap frequency as a function of the inter-well potential height. Experimental data (circles) 
are compared to the values calculated with Eq. (grey triangles) and to a numerical 
simulation of the ID GPE (open triangles). 



in the lattice and to the frequency uj x of the condensate oscillations in absence of 
the periodic field by the relation 

AhN T ( u x ' 



ic = ^r — • (16) 

Figure (right) shows the experimental values of the oscillation frequencies 
together with the result of a variational calculation. It must be noted that, due to 
mean field interactions, in our system only for potentials higher than ~ Er a bound 
state exists in the lattice. 

3.6 Expansion inside a moving ID optical lattice 

We now discuss a different experiment wh ere we study the e xpansion of a con- 
densate inside a moving ID optical lattice (jFallani et aZl 12003^ . This experiment 



allows us to load the condensate with different quasi-momenta q in the periodic 
structure realized by the optical lattice. We will show how it is possible to adopt a 
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completely different language, namely that of band structures in solids, to describe 
the condensate motion. This is possible because the momentum spread in the con- 
densate is so small that the behavior of the entire atomic cloud can be explained 
as that of a single particle. From solid state physics it is well known that in the 
presence of an infinite periodic potential the ene rgy spectrum of the free particle 
is modified and a band structure arises (see e.g. Jones and March] . Il986h . In the 
rest frame of the lattice the eigenenergies of the system are E n (q), where q is the 
quasi-momentum and n the band index. According to band theory, the velocity in 
the n-th band is v n = Tr x dE n ldq and the effective mass is m* = h 2 (d 2 E n /dq 2 )~ 1 . 
The effective mass can be negative for a range of quasi-momentum and this has 
been recently recogni zed as a possibility of realizing bright soliton s in BEC with 
repulsive interactions ( Lenz et al . 19941 : Konotop and SalernoL 20021 : HilligsO et al. 



2002; lEiermann e~l 



2003h . 



In this experiment we first produce the condensate in a pure harmonic trap, then 
we switch off the magnetic harmonic potential let the BEC expand for 1 ms and we 
switch on a moving periodic potential. After 1 ms of expansion the density of the 
condensate decreases enough to neglect the non linear term in the Gross-Pitaevskii 
equation describing an interacting BEC. This means that, as a first approximation, 
we are allowed to consider the BEC as a linear probe of the periodic potential 
energy spectrum. The moving periodic potential is created by the interference of 
two counterpropagating laser beams with a slightly different frequency and blue 
detuned 0.5 nm from the D2 resonance at 780 nm. The two beams are obtained by 
the same laser and are controlled by two independent acousto-optic modulators. The 
resulting light field is a standing wave moving in the laboratory frame with a velocity 
vl = XAu/2, where Av is the frequency difference between the two laser beams. 
In our experiment, we typically vary the optical lattice velocity between and 2vr 
where vr = hki/m is the recoil velocity of an atom absorbing one lattice photon 
and corresponds, in the frame of the band theory, to the limit of the Brillouin zone. 
We switch on the moving optical lattice adiabatically by ramping the intensity of 
the two laser beams in 2 ms. This ensures we are lo ading the condensate in a Bloch 
state of well-defined energy and quasi momentum ([Denschlag et al. . 2002). We let 



the condensate expand in the lattice and after a total expansion time of 13 ms we 
take an absorption image of the cloud along the radial horizontal direction looking 
at the position and dimensions of the condensate inside the optical lattice. From 
the position after the expansion, we extract the velocity of the condensate inside the 
optical lattice. In particular, we repeat the experiment for different velocities of the 
lattice and compare the position of the expanded condensate inside the lattice with 
the position of the condensate expanded without the optical lattice. Let us call Az 
the difference in position along the axial direction. Then, the velocity of the BEC 
inside the optical lattice is given by v = Az/ At — vl, where At is the expansion 
time inside the lattice. 

In Fig. |H1 we show the results obtained for the velocity of the condensate as a 
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Figure 8: Velocity of the condensate inside the optical lattice as a function of the 
quasi- momentum q in units of the recoil momentum q^ = hk^. The open circles cor- 
responds to data obtained with V op t = 1.3Er and the filled circles to data obtained 
with V op t = 3-8Er. The dashed and dotted lines are the correspondent curves given 
by the band theory. 



function of the quasi-momentum q for two different values of the lattice potential 
depth. The experimental data points are compared with the theoretical results 
obtained from the band theory and show a very good agreement. With an adequate 
sampling of the velocity we can extract the effective mass values given by dv/dq. 
The results for an optical potential depth of 1.3 Er are shown in Fig. As we 
increase the lattice velocity (corresponding to increasing the quasi-momentum q) 
the effective mass rapidly increases and between 0.7 q and 0.8 q it first becomes 
infinite positive and then negative. 

The consequence of the strong variation of the effective mass is expected to 
consistently modify the expans ion of the condensate along the axial direction 
( Massignan and Modugnd l2003h . As a matter of fact, the effective mass enters 
the diffusive (kinetic) term in the Gross-Pitaevskii equation. 

In Fig. EI] we report the radii of the condensate measured as a function of the 
quasi-momentum after the expansion inside the optical la ttice compared to numeri- 
cal p redictions based on an effective ID theoretical model ( Massignan and ModugnoL 
The axial radius (filled circles in Fig. I10|) decreases until the quasi-momentum 
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Figure 9: Effective mass of a condensate moving in an optical lattice of 1.3Er as 
a function of the quasi-momentum. The data points correspond to the values ex- 
tracted from the measured velocity and the solid line is the corresponding theoretical 
prediction of the band theory. 



reaches qR. This is first due to the increase of the effective mass (causing a slower 
expansion) and then by the fact the the effective mass becomes negative (causing a 
contraction of the axial direction during the time spent in the optical lattice). 

When q > qR, the effective mass becomes positive again but with a value smaller 
then the real mass m. As a consequence the expansion becomes faster in this 
region of quasi-momenta. In Fig. we also report the measured values of the 
radial dimension of the BEC. A deviation from the expansion without optical lattice 
(dashed line) is observed also in this direction for q < qR, even if this dimension 
is not directly affected by the presence of the lattice. This is consistent with the 
theory (continuous line) and can by explained in terms of a coupling between the 
axial and the radial dynamics. 

For q < qR, the compression along the lattice direction increases the mean- 
field energy and causes a faster radial expansion. Instead, when the condensate 
is loaded with q > qR, the axial expansion is enhanced (0 < m* < m) and the 
residual mean- field energy is further reduced, causing a suppression of the non-linear 
coupling between the axial and radial dynamics. This behaviour is evident in the 
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Figure 10: Axial and radial dimensions of the condensate after the expansion in 
an optical lattice V pt = 2.9-Er. The experimental points (filled and open circles) 
show the Thomas-Fermi radii of the cloud extracted from a 2D fit of the density 
distribution. The dotted lines show the dimensions of the expanded condensate in 
the absence of the optical lattice. The continuous and dashed lines are theoretical 
calculations obtained from the ID effective model. 



absorption images reported in Fig. 1111 where we show the shape of the condensate 
expanded without optical lattice (image a)), and with an optical lattice of 2.9 En 
and respectively quasi-momenta q < qn (b)) and q > qu (c)). In the first case, a 
contraction along the axial direction is accompanied by a faster expansion along the 
radial direction, while in the second case the condensate expands faster in the axial 
direction. 



3.7 Experiments beyond the GP equation 

More recently experiments are starting to explore the possibility of cre ating and ma- 



nipulat ing pure quantum states of many atoms. In a beautiful work, iGreiner et al 



2002a) haw has loaded a BEG in a three-dimensional optical lattice. By adiabat- 



ically increasing the periodic potential, it was possible to enter in a regime where 
the tunnelling between adjacent wells was comparable with the atomic interaction. 
In this situation, quantum fluctuations of the atom number in each well become 
relevant and the system can no longer be described by the GP equation. The sys- 
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a) 



b) 



c) 



Figure 11: Absorption images of the expanded condensate. From left to right: 
(a) normal expansion of the condensate without lattice; (b) axial compression in a 
lattice of 2.9 and vl = 0.9vr; (c) enhanced axial expansion in a lattice of 2.9 Er 
and vl = 1-Ivr. 



tern undergoes a Mott-insul ator quantum phase transition to a s tate where the atom 
number in each well is fixed (|.Taksch et al!ll998l : lslchdevl .ll999). In such a situation 
no interferogram can be recorded when the atoms are released from the potential. In 
addition, it is not possible to move the atoms through the array any longer until the 
potential energy difference between two adjacent wells equals the energy necessary 
to add an atom to an already occupied site. 

This experiment, aside from showing the first clear quantum phase transition, 
opened the path for a new set of experiments where it was possible to show col- 
lapses and revivals of the matter wave field , as non-trivial quantum states were 



formed in the lattice Bremer et all B . Later, in the same system, it was 



possible to create massive entanglemen t between distant a toms in the lattice by 
the coherent manipulation of collisions ( Mandel et al. L 20031) and, more r ecently, a 



Tonks-Girardeau gas was observed for the first time (Parcdcs et an . 12004) 



4 Limitations of the GP equations from first-principles 
arguments 

The purpose of this short section is to expose some limitations of the GP equations 
by appealing to fi rst-principles th eory. The most fundamental approach seems to us 
to be that due to LeggettJ (|2003h . The arguments outlined below therefore largely 
follow his article, though because of space limitations our discussion will be more 
qualitative than his original study. 

His starting point is what he terms the formulation of a 'pseudo-paradox' in 
the theory of a dilute Bose gas with repulsive interactions. His 'paradox' can be 
stated as follows. The usual ground-state energy in the GP approximation set out 
in Section [2] above is lower than that in the Bogoliubov theory. 1 Thus, by standard 
variational arguments, the GP answer should be an improved approximation to 



1 Concerning the calculation of the ground-state energy of a homo geneous Bose-Einste i n con - 
densate beyond Bogoliub ov theory, we refer the reader to the works of lLieb and Yngvasonl 
IWeiss and Eckardd i20o3) . and references given therein. 
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the Bogol iubov re s ult, w hich contradicts the 'established wisdom' concerning this 
problem. LeggettJ ( 20031 ) then resolves the 'paradox' by a correct transcription of 
the two-body scattering treatment to the many-body problem. He stresses that this 
resolution has not to do with spurious ultraviolet divergences resulting from the 
replace ment of th e true interatomic potential by a model (^-function pseudopotential. 
Instead, iLeggetti points out that it comes from an infrared divergence which has as 
a consequence, first of all, that the GP app roximati on actually has no well-defined 
many-body wave function underpinning it. Leggettl goes on then to show that the 
'best' attempt to construct an approximate variational wave function always results 
in a ground-state energy which either exceeds, or at the very best is equal to, that 
given by the Bogoliubov approximation. 

This all prompts us to follow this brief section with an account of a proposal 
bv lAngilella et ol\ ((20041 ) to replace the static GP differential Eq. Q by an integral 
equation (SecEI). The starting point of their derivation is the Bogoliubov-de Gennes 
theory, and therefore some physical background to their theory will first be given in 
Section immediately below. 



5 Bogoliubov-de Gennes equations: physical back- 
ground 

A key ingredient of the GP equation is a spatially varying order parameter ^(r), 
describing the inherently inhomogeneous density profile of a BEC condensate. Such 
an assumption is usually not necessary for a pure, conventional (i.e., s-wave) super- 
cond uctor, for which BCS theory applies with a constant order parameter A (see, 
e.g., 



Schrieffer, This is not the case for inhomogeneous superconductors, such 



as alloys or metals in the presence of impurities, both diluted and isolated, where 
a position dependent order parameter A(r) is required. Such a generalization is 
also relevant to describe boundary effects in superconductors, or the effect of a non- 
uniform magnetic field, as is given by a position dependent vector potential A(r). 
The set of coupled equations which relate the inhomogeneous pairing potential A(r) 
and the self-consistent potential U(r) experienced by single particles in the super- 
conductor can be derived within a mean-fi eld approach originally due to Bogoliubov 
( Bogoliuhov etail ITflRflt IHb Oermesl Rflfi^. 



One usually starts with the effective, mean-field Hamiltonian 

H = y"dr$,(r)[Wo + tf(r)]$ Q (r) + J dr [A(r)^}(r)^{(r) + H.c] (17) 

generalizing the BCS effective Hamiltonian to the inhomogeneous case. Here, ^a(r) 
[^>Q.(r)] is a creation [annihilation] Fermion field operator (with understood sum- 
mations over repeated spin indices), A(r) is a spatially varying pairing potential, 
self-consistently taking into account for the electron-electron interaction, which is 
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here assumed to be described by a spin independent, strongly local potential, as in 
Eq. @. In Eq. p7|). the kinetic term reads 

w »=--£( v -s a ) 2 -"' < is > 

with [i again the chemical potential and A(r) the electromagnetic vector potential, 
and U(r) is the self-consistent potential experienced by single particles. 

One may conveniently expand the field operators over a complete basis set as 



tp a (r) = ^2 0„(r)c„ a , 



(19) 



where c va , ct a obey the usual Fermion anticommutation rules, and <f> v (r) are eigen- 
functions of H.q with eigenvalues £ v (measured from the Fermi level). Depending on 
the particular symmetry of the proble m under study, 4>,Ar) may be plane waves, or 
more complicated functions (see, e.g., iGvgi and Schliiterl . Il99ll . for an instance of 
axially symmetric systems). 

Since the effective Hamiltonian, Eq. (|17[). is quadratic in the field operators, it 
can be diagonal ized by means of a Bogoliubov-Valiatin canonical transformation 
( SchrieffeiJ . ll96l l. )ut now allowing for spatially varying coherent factors u u (r) and 
v v (r): 



^ T (r) 



V 



(20a) 
(20b) 



with jt a pjw] Fermion creation [annihilation] operators for the quasiparticles in 
the superconducting state. 

One then requires that Eqs. (12011 diagonalize the effective Hamiltonian, Eq. (|17j) . 

as 

H = E +J2wta%a, (21) 
va 

where Eq is the energy of the superconducting ground state. One then easily finds 
i lua\ — — ti/'Jva and [^,7^0.] — tv^fta- On the other hand, from Eq. (|17|) . one has 



-Hih(r) - A(r)$(r) 
■«^(r)+A(r)&r) 



(22a) 
(22b) 



where TL = TCo + U(r), and making use of the transformations Eqs. (120H one even- 
tually obtains 



e u u v (r) = Hu u (v) + A(r)«„(r) 
e^(r) = -H\(r) + A*(rK(r), 



(23a) 
(23b) 
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whi c h are the static Bogoliubov-de Gennes (BdG) equations ([Bogoliubov et al. 
19591 : Ide Gennesl . Il966h . 



Eqs. (|23j) can be written in compact matrix form as an eigenvalue problem for 
the coherence factors: 



H A(r)\ (u u (r)\ _ (u v {v] 
A*(r) -H*J \v v (r)J ev \v u {v) 



(24) 



Within such a mean- field approach, the pairing potential A(r) and the single- 
particle potential U(r) are identified with the self-consistent statistical averages 

A(r) = g(^(br)^(br)) (25a) 
U(r) = g(^\(br)^(br)), (25b) 

which, on making use of Eqs. (|20|). become: 

A(r) = -gJ2<(r)u v (r){l-2f v ) (26a) 

V 

U(r) = gJ2[K(r)\ 2 fu + \vAr)\ 2 (l-U)], (26b) 

V 

where f u = [exp(e u /kBT) + is the Fermi function evaluated at e v . Thus, su- 
perconductivity in an inhomegeneous system is governed by the eigenvalue problem 
given by the BdG Eqs. I)23j) . and the two self-consistency conditions, Eqs. (|26|) . 
In analogy with Eq. 0, also the BdG equations can be generalized to the time 
dependent case. 

One possible variational derivation of the BdG equations r ests on the generaliza - 
tion of density functional theory (DFT) to superconductors ( Oliveira et al. . 1988^ . 
There, Eqs. (|23jl can be identified with the Kohn-Sham equations for an appropriate 
Hohenberg-Kohn functional of both particle density 

n(r) = J>+(r)4(r)> (27) 

(J 

and the superconducting order parameter 

x(r,r') = (^ T (r)^(r')), (28) 

here generalized to take into account nonlocal effects. Such an approach has been 
recently employed to evaluate the superconducting condensation energy of the ho- 
mogeneous electron ga s with anisotropic pairing potent ials, with angular momentum 
t ranging from 1 to 9 (j Wierzbowska and Kroghl l2005h . 
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6 Integral equation transcending static GP equation 
taking Bogoliubov-de Gennes theory as starting point 



In a very stimulating recent contribution, Pieri and Strinatl (|2003h (referred to as PS 



below) have 'derived' the non-linear GP differential equation for condensed bosons 
by taking as their starting point the Bogoliubov-de Gennes equation for superfluid 
fermions. 

The purpose of this section (see also Angilella et al. , 2004^ is to demonstrate that 



one can generalize the zero-temperature differential GP equation while remaining 
within the original framework of PS, an integral equation formulation then resulting. 
The framework of PS is provided by the coupled integral equations involving Green 
functions G21, G\\ and G . The equations are: 

Gn(r,r'; Ws ) = G D (r, r'; u.) + jdv"G {v, r"; ^)A(r")G 21 (r", r'; u s ), (29a) 

G 2 i(r,v';u s ) = - Jdv"G {v" ,v; -u s )A*(r")G n (r" , r';u s ), (29b) 

where uj s = {2s + 1)tt/P (s is an integer) is a fermionic Matsubara frequency, 
(3 = l/k-sT, Gn is the normal and G21 is the anomalous single-particle Green 
function. The third Green function appearing in Eqs. I|29j). namely Go, satisfies the 
equation 

[iu s - H(r)]G (r, r'; u s ) = 6(r - r'), (30) 
where the single-particle Hamiltonian H (r) is defined by: 

H(r) = -^-V 2 + V(r)-^ (31) 

\i being the Fermionic chemical potential. As PS stress, Eqs. 1)29(1 . when taken 
together with the self-consistency equation for the gap function: 

A*(r) = ^^G 21 (r,r;^), (32) 

s 

are entirely equivalent to the Bogoliubov-de Gennes equations that describe the 
behavior of superfluid fermions in the presence of an external potential. Equa- 
tions ((29H32() define what we have termed the original framework of the PS study. 
The constant Vo < entering Eq. 1(32(1 arises from the contact potential Vo<5(r — r') 
assumed by PS to act between fermions with opposite spins. We also retain here 
their use of the ratio A(r)/// as an expansion parameter which allows the rapid trun- 
cation of such series, which then leads for strong coupling to an integral equation 
for the gap function 

-^-A*(r) = J driQ(r,ri)A*(n) 

+ J dr 1 dr 2 dr 3 -R(r,r lj r 2 ,r3)A*(r 1 )A(r 2 )A*(r 3 ), (33) 
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where R is written explicitly in terms of G Q (r, r\;uj s ) in Eq. (15) of PS. However, 
as will emerge below, it is the non-local kernel Q(r,r') which is at the heart of the 
present study. In terms of the Green function G entering Eq. (|3U|) . Q(r,r') is given 
by [PS: Eq. (14)]: 

1 

P 



Q(r,r') = -^G (r',r;-a; s )G (r , ,r; Ws ). (34) 



We take the integral equation (|33|) for the gap function as the starting point of this 
section. For our purposes below, it is then crucial to gain insight into the kernel 
Q in Eq. (|34[). and in particular to carry out the summation explicitly over the 
Matsubara frequencies oj s . 

To gain orientation, let us first perform this summation when the external po- 
tential V(r) is set to zero in Eq. 1)30(1 . Having achieved this summation, we shall 
present a general m ethod to allow the su m over uj s to be achieved for V(r) ^ 0, 
using earlier work of Stoddart et ol\ (Ej 



Returning to the explicit form of Q(r, ri) given in Eq. (|34[) above, it is natural 
to study first the translational invariant, free-electron limit of Eq. 1)34(1. say Qo(t), 
with r = |r — ri|, which is obtained by 'switching off' the one-body potential V(r). 
This amounts to replacing G in Eq. 1)34(1 with the free-electron Green function G . 
For the Fourier transform of Qo(r), we formally find 

4 _ f * i-Mf*-*)-M<x) (35) 

J (2vr) d Ck-k' + 4k' 

where £k = k 2 /2m — [i and np(£) is the Fermi-Dirac distribution function. However, 
it should be noted that, in three dimensions, Eq. 1)35(1 contains a divergent contri- 
bution at large wave-numbers, which implies a divergent behavi or of Q n (r) at sma ll 
distances r. Indeed, we find the asymptotic expansion (see also Alexandrov . 2003h : 



where r' = kpr, kp is the Fermi wave-number, defined by fj, = kp/2m, f3' = fifj,, and 
a = r'n/ (3' . Fig. 1121 shows then our numerical results for r' 4 Q function of 

r', for several temperatures {$' = 1 — 30). 

Following Stoddart et al. ( 1968^ . the canonical density matrix C(r, r',/3) is de- 
fined by 

C(r,r',P) = J2MrM(r')e- l3e \ (37) 



wher e (5 = 1/ksT. Within the perturbative approach of March and Murravl l|l96Cll . 



196lh . with plane waves as the unperturbed solution, the canonical density matrix 



can then be written to all orders in the external potential V(r) in terms of the 
free-particle canonical density matrix given by 

C (z, (3) = (2ttP)- 3 / 2 exp(-z 2 /2P), (38) 



30 



G. G. N. Angilella, S. Bartalini, F. S. Cataliotti et al. 




Figure 12: Solid lines show r' 4 Q (r'), where Q(r) is defined by Eq. ()34[). as a function 
of r' = k-pr, for several temperatures, given by (3' = = 10 — 30 (bottom to 
top). Dashed lines are the asymptotic expansion Eq. (|36|) . Units are such that 
fcf/(4£t) = 1. Redrawn from lAngilella et q/J(|2004h . 



MS 



C(r,T O ,0) 



dzzCo(z,/3)f(z,r,T ), 



(39) 



where / satisfies the integral equation ( Stoddart et q/l ll968): 

/( Z) r,r ) = i^-|r-r |)- / dr t ^ f(z - |r - r x |, n, r ). (40) 
z J Z7r|r — n| 

The des ired Green function G n is then to be obtained from / entering Eqs. (|39|) and 
(liO]) as (jStoddart pi ll96Sh 



G (r,ri; k) 



dz zG (z;k)f(z,r,ri), 



where 



G (z;k) 



Akz 



4-7TZ 



(41) 



(42) 



One may also take advantage of the expression in Eq. (|41|) of G in terms of G Q to 
rewrite the kernel Q(r,ri) defined by Eq. (|34|) as 



/■oo 

Q(r,ri)= / dzidz 2 z 1 Z2f(zi,ri,r)f(z2,ri,r)Q (zi,z 2 ) 
Jo 



(43) 
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where the Fourier transform of Q (zi, z<i) is given by 

Qo(k 1 ,k 2 ) = 1 - n f^; nF( ^ ) . (44) 

Ski + ?k 2 

Hence, the sum over Matsubara frequencies has still been carried out in the presence 
of an external potential V(r) entering Eq. (|4Uj) for the function /. 

Because of current interest in harmonic confinement in magnetic traps at low 
temperatures, let us illustrate the rather formal Eqs. (|39|) and (|4U[) when the ex- 
ternal potential V(r) has the explicit isotropic harmonic oscillator form in three 
dimensions, namely 

V(r) = -muV. (45) 

Following the pioneering work of ISondheimer and Wilsonl (jl95ll ) on free electrons 
in a magnetic field, the diagonal element C (r,r,0) when V(r) is given by Eq . (|45() 
takes the form (see e.g. March et all 19951 . p. 27; see also Howard et a,L . 2003) 



C ^ = (ssT (ffl) 3 ' 2 ^ (-^tanhl^) , (40) 

which is the so-called Slater sum of quantum chemistry (Fig. 113(1 . 

From Eqs. (|38f) and (J39|). performing the substitution t = z 2 /2, it then follows 
that f(z,r, ro) can be expressed as the inverse Laplace transform 

/(*, r, r ) = (2^) 3 / 2 £- 1 [ S -^ 2 C(r, r , s' 1 )] (47) 

where (t, s) are conjugate variables with respect to the Laplace transform. 
Within the Thomas-Fermi (TF) approximation, we take: 

C TF (r,r,/3) = 1 exphm*)], (48) 

which is plotted also in Fig. ^3 f° r ^( r ) given by Eq. (|45jl. For the value of 
shown, the TF form Eq. (|48() is seen to be a useful approximation to the exact 
result, Eq. (@BJ). Inserting Eq. (gHJ into Eq. (gTJ) we find 



/TF(,,r,r) = S M - fflj!^], (49) 

where Ji denotes the Bessel function of the first kind and order one. 

Fig. El shows f(z, r, r) as a function of r for fixed z, as obtained by numerically 
performing the inverse Laplace transform in Eq. (|47|) for the harmonic potential 
case. The regular contribution to the analytic result for the Thomas-Fermi approx- 
imation, Eq. (|49|). is also plotted for comparison. The similarity in shape between 
approximate and (numerically) exact results for this harmonic confinement model 
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Figure 13: Shows diagonal element of the canonical density matrix C(r, r, /?), 
Eq. (|46|). and its Thomas- Fermi approximation, Eq. (|48|). as a function of r, for 
/3 = 0.5. En ergies are in un i ts of faj, while lengths are in units of (h/muj) 1 ^ 2 . 
Redrawn from Angilella et al. ( 20041 ) . 



seems to us rather remarkable. After this model test of a TF-like approximation 
invoked by PS, we return to the general case, based on the exact result Eq. (|4"3"|) for 
the kernel Q(r, r±). 

Then, we invert the argument of PS but still use a further essential assumption 
of their study, namely that the condensate wave function ^(r) entering the Gross- 
Pitaevskii equation is related to the gap function A(r) by 

Here, in the strong coupling limit, and following PS, ~ (2m\n\)~ l l 2 represents 
the characteristic length scale for the non- interacting Green function G , equal to 
G above when V(r) is put equal to zero. 

Given the validity of this PS assumption, Eq. (|5U|). we then rewrite Eq. ()33|) as 
an equation for <I>(r): 

-^r**to = / <friQ(r,n)**(n) 

+ ^2 / dr 1 dr 2 dr 3 R(r,r u r2,r 3 )^*(r 1 )<P(r 2 W(r 3 ). (51) 
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Figure 14: Shows diagonal f(z, r, r) corresponding to the harmonic potential, as 
given by the inverse Laplace transform, Eq. (|47|). as well as the regular part of its 
Thomas-Fermi approximatio n, Eq . (14911 . as a funct ion of r, for fixed z = 0.5. Units 
as in Fig. 1131 Redrawn from Angilel la et al. ( 20041 s ) . 



This then is the proposed generalization of the Gross-Pitaevskii equation, but with 
Q(r,ri ) to be calculated more accurately than by the Thomas- Fermi-like assump- 
tion of iPieri and Strinatil ((20031 ) . via Eqs. (j33|) and @Qj) . 



While Eq. (|51|) is a direct consequence of the above arguments, it remains an 
expansion in <£, in suitable reduced form. Therefore, a first attempt to simplify this 
Eq. (|51|) is to retain the approximation given by the Pieri-Strinati approach in the 
'smallest' term involving 0(<I> 3 ) on the right-hand side of the basic Eq. (|51|) . Thus 
one reaches the (still non-local) equation for the condensate wave function <I>(r): 

-^-$*(r) = f dn Q(r, n)*( ri ) - ^ma||$(r)| 2 $(r). (52) 

For sufficiently small spatial variations in the condensate wave function ^(r) in 
Eq. (|52j) . the basic nonlocality can be removed by Taylor expanding ^(ri) around 
the position r in the integral term. This then characterizes the problem in terms 
of 'partial moments' of the kernel Q(r,ri), namely / Q(r, ri)dri and / Q(r, ri)|r — 
ri| 2 dri. Such partial moments then enter the original GP equation, as stressed by 
PS. 

In summary, we propose the retention of the non-local kernel Q(r,ri) as in 
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Eq. H51|) above, since the sum over Matsubara frequencies in Eq. (|34|) has been 
performed in Eq. Q43|) . which is a central result of the present study. However, 
in the terms of 0( < I )3 ) in Eq. (|51|). a sensible starting point is to follow the PS 
approximation displayed in Eq. (|52[). 

As to future directions, evaluation of the non-local kernel in Eq. (|43j) for other 
external potentials than the harmonic case in Eq. (|45[) is of obvious interest. For 
this latter model, though our Fig. 1141 conside rs the diagonal eleme nt of f(z, r, ri), 
the off-diagonal form of C(r,ri,/3) is known ( Howard et al\ . 2003h . and numerical 
Laplace inversion to obtain /(z, r, ri) is entirely feasible. Then Q(r, ri) can be 
obtained, though of course numerically. 

The GP equation is valid in the strong-coupling limit of superfluidity. It has to 
be stressed that in the weak coupling limit one can also derive a Ginzburg-Landau 
equation starting from the Bogoliubov-de Gennes equations. We note specifically 
in this context that the derivation of the Ginzburg-L andau equation in th e weak - 
coupling limit for the harmonic trap was presented by Baranov and Petrov ( 1998h . 
The results presented in this section are also relevant to the weak-coupling limit of 

superfluidity. 

Finally, we return to the foundations of the Gross-Pitaevskii e quation b y Leggettl 

which we summarized very briefly in Section above. iLeggett concludes 
that there is no correlated many-body wave-function underlying their original equa- 
tion. It will be interesting for the future to know whether the non-local versions of 
Eqs. (|5l"|) and (|52|) proposed here are still subject to this limitation. 



7 Summary and suggestions for further study 

In summary, the GP Eqs. (HJ and © have been considered in relation to a va- 
riety of experiments on inhomogeneous condensed Bosons in Sections |2] first and, 
predominantly, in Section [3 

Sections BHB1 are then theor etical in emphasis. The short Section 0] summarizes 



the arguments of lLeggett ( 20031 ) , which reveal that there is no underlying many-body 



wave function for the static GP Eq. ©. Therefore, instead of the GP differentia l 
equation, SectionsEJandElpresent an integral equation theory ( Angilella et al . 2004h 



for which the starting point is the Bogoliubov-de Gennes theory. 



7.1 Suggestions for further study 

Though a non-local theory given in Section H3 in the form of an integral equation, 
trascends the static GP differential Eq. @, no proof has, as yet, been given that this 
integral equation theory has, underlying it, a many-body wave function. It seems 
that this, therefore, in view of Leggett's criticism of the foundations of the static 
GP Eq. is an area of considerable interest for further theoretical work. 
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O'Dell et flZJ (l2004h have developed, within the Thomas-Fermi regime, the hy- 



drodynamics of a trapped dipolar BEC. As they stress, such a BEC, whose par- 
ticles interact via dipole-dipole coupling, constitute an example of a superfiuid 
with long-range and anisotropic interparticle forces. This contrasts with the cus- 
tomary BECs which are characterized by isotropic interacti ons with range muc h 



smaller than the average interparticle interaction (see e.g. Dalfovo et al. . 1999) 



While alkali atoms have small dipole-dipole interactions, O'Dell et al\ suggest that 
chromium is a promising case for a dipolar BEC, due to its large magnetic moment 
of 6 Bohr magnetons, and there has been p rogre ss in cooling i t towards degeneracy 
( Schmidt et all 12003; H ensler et all I2003T ). As lO'Dell et al\ also note, molecules 



can possess larg e dipole moments. They suggest that advances in the cooling of 
polar molecules dWeinstein et all llQQfit iBethlem et oil. Il999l Eoool l2002h . photoas- 
sociati on of ultracold heteroriuclear molecu l es dMancini et all l2o7)4r) . an d molecular 
BECs dCreiner et all 1200.4 I.Tochim et all 1200.4 IZwierlein et all 12003^ may soon 
lead to superfluids where dipolar effects play a major role. 

O'Dell et al\ (2004) consider in their Letter a harmonically trapped BEC with 



dipole-dipole interactions as well as short-range s-wave scattering. They note that 
in the Thomas-Fermi limit, where the zero-point kinetic energy of the atoms in 
the trap can be neglected compared with the interparticle interaction energy and 
the trapping potential, the collective dynamics of a BEC may b e treated by th e 
collisionless hydrodynamic theory of Bose superfluids at T = ( Stringarl . 1996h . 
The effect of the dipolar interactions is to introduce non-locality into the already 
nonlinear hydrodynamic equations. 

In rel ation to the focus of the present review on the GP equations, it is to be em- 
phasized ( O'Dell et al. , 2004) that the dominant interactions in the ultracold gases 
available at the time of writing are asymptotically of van der Waals character, with 
an r -6 decay, and are short-ranged, compared with the mean interactomic distance. 
As we have seen above, within the mean-field regime of the GP equations these 
interactions are represented by the pseudopotential model gS(r) = (Aira s H 2 /m)6(r) 
[cf. Eqs. © and ©]• By employing a Fano-Feshbach resonance, it is possible to 
adjus t the value of a, to yield repulsive (positive) and attractive (negative) val- 
ues |lnouve et al. , 199£lh . Appealing to the analogy of nuclear magnetic resonance 
techniques, di pole-dipole interactions can also be controlled by rapidly rotating an 
external field (|Giovanazzi et al. . 2002h . The interactions vanish when the rotation 



corresponds to the so-called magic angle. 



To summarize here the main points made bv lO'Dell et all (|2004l ). the versatility 



of quantum gases is promising for the study of the role of interactions in superflu- 
idity. They present solutions of the dipolar superfiuid hydrodynamic equations in a 
harmonic trap: the condensate density is parabolic as in the pure s-wave problem 
but now with modified radii. 

Finally, and returning to the Fano-Feshbach resonances already referred to 
( Fanol . 193.4 Feshbacll . 19581 ) . there remains the important question for theory of the 
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crossover of Bosonic to Fermionic super fluidity (see, e.g., Randerial. 19951) . following, 
for example, the pioneeri ng; studies of Nozieres and Schmitt-Rinkl (jl985T ) (see also 



Timmermans et all 12001). Note that the cold-atom technology that has resulted in 



the observation of dilute-gas BECs enables the creation and subsequent studies of 
physical properties of novel superfluids having a degree of flexibility adding to the 
already powerful methods of traditional low temperature physics. This therefore 
seems a further, and very attractive, direction for further studies. 
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